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DECREASING  DG  AND  DIRECTION  DISPERSION 
OP  MAIN  MAXIMUM  OF  BEAM  PATTERN  OP 
SECTIONAL  TRAVELING  WAVE  ANTENNA 


V.  I.  Zamyatin  and  L.  G.  Korniyenko 

Summary 


Formulas  are  found  and  analyzed  to  reduce  the  directive  gain 
and  direction  dispersion  of  the  main  maximum  of  the  beam  pattern  of 
a  sectional  traveling  wave  antenna  when  there  are  errors  of  the  wave 
number  in  the  system  and  phase  errors  at  the  sites  of  contact  between 
the  sections.  A  condition  of  ootimality  of  sectioning  and  a  con¬ 
dition  m  which  the  limiting  directive  gain  is  missing  are  obtained. 
The  effectiveness  of  the  sampling  of  phase  errors  is  studied. 
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Introduction 

Production  conditions  force  us  to  assemble  large  antennas  made 
of  individual  sections. 

Since  the  dimensions  of  the  sections  are  generally  small ,  there 
is  a  possibility,  first  of  all,  of  manufacturing  them  fairly  accurately, 
and  secondly,  of  setting  up  comparatively  simply  the  measurement  of 
the  amplitude-phase  distribution  (APD)  in  .the  aperture  of  the  sections 
and  its  deviation  from  the  calculated  value  (errors). 

The  latter  circumstance  permits  a  more  rational  approach  to  assembly 
of  the  antenna,  that  is,  use  of  the  methods  of  assembly  (sectional)  which 
reduce  the  effect  of  random  errors  in  APD  on  the  antenna  parameters. 

A  number  of  works  have  covered  a  statistical  analysis  of  sectional 
antennas  [1-3]. 

Publication  [1]  has  studied  the  statistics  of  the  field  of  inphase 
sectional  antennas.  Publications  [2-3]  have  investigated  the  charac¬ 
teristics  of  linear  [2]  and  two-dimensional  [3]  sectional  traveling 
wave  antennas  (STWA)  with  different  methods  of  sectioning. 

Publication  [21  in  particular,  has  exatained  the  STWA  characteristics 
with  random  and  selective  methods  of  assembly.  It  has  introduced  an 
important  parameter  which  characterizes  the  economic  efficiency  of 
production.  Publication  [3]  has  found  the  average  beam  pattern  of  the 
system  which  consists  of  a  set  of  linear  STWA  with  random  and  one-sided 
method  of  assembly,  with  sampling  of  errors  at  the  end  of  each  section. 

One  of  the  main  results  of  these  works  is  the  conclusion  that  with 
an  increase  in  the  number  of  sections  ,  the  distortions  of  the  beam 
pattern  (BP)  in  antenna  N  diminish.  This  conclusion  is  correct  in  the 
case  where  the  junction  sites  of  the  sections  are  ideal.  In  real 
antennas ,  at  the  junction  sites ,  inhomogeneities  naturally  arise  which 
result  in  random  errors  in  the  APD.  In  this  case,  enlargement  of  N 
results,  on  the  one  hand,  in  a  decrease  in  errors  governed  by  fluctuations 
in  the  wave  number  ,  and  on  the  other  hand ,  in  an  increase  in  the  errors 
associated  with  the  nonideal  nature  of  the  junctions.  The  same  situation 
develops  when  the  errors  at  the  end  of  each  section  are  not  sampled 
completely,  but  with  a  certain  error.  ^ 

In  these  cases ,  a  certain  optimal  solution  evidently  exists  in 
selecting  the  number  of  sections. 
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In  addition,  the  sanding  of  phase  errors  to  zero  at  the  end  of 
each  section  which  is  examined  in  [3]  can  be  difficult.  It  is  therefore 
important  to  study  a  more  general  case  where  the  error  sampling  is  made 
only  through  M  sections. 

This  work  studies  the  statistics  of  STWA  with  the  same  methods  of 
sectioning  as  in  [2],  but  first  of  all,  with  sampling  of  the  phase 
errors  to  zero  through  M  sections  ,  and  secondly ,  with  consideration  for 
the  phase  errors  which  develop  at  the  sites  of  the  section  junctions. 

For  convenience,  the  same  letter  designations  are  adopted  as  in  [2]. 


Initial  Correlations 

Assume  that  a  STWA  of  length  L  is  formed  of  N  sections  equal  to 
length  1  =  taken  from  the  original  group  of  sections. 

When  the  antenna  has  phase  errors ,  separate  realization  of  the 
beam  pattern  for  output  looks  like: 


F(u,  a)  = 


e 


i  ^  (jii-i) 


P  <*> 

e  *  dy  , 


(1) 


where  An(y) --amplitude  distribution  in  n-th  section  of  STWA, 

0  (y)  =  ^SiM  --fluctuations  in  the  phase  in  the  n-th  section,  classified 

— -  0  as  a  certain  quadratic  value  o , 

u=  ^(kg-k  cos  0) --generalized  angle, 

k,  kQ--wave  numbers  in  system  and  in  free  space  respectively. 

For  simplicity,  we  will  consider  that  the  amplitude  distribution 
in  the  antenna  aperture  is  uniform,  and  the  phase  velocity  of  the  wave 
in  the  system  is  greater  than  or  equal  to  the  speed  of  light.  Then, 
in  the  absence  of  errors ,  the  ^irection  of  the  main  maximum  Uq=0  and 
the  beam  pattern  F(u,0)=L^  — •  When  there  are  small  phase  errors, 
the  directional  drift  of  the  main  maximum  Au  is  determined  by  the 
ratio  (A- 3)  with  b=o  and  r.4««  and  equals*): 


f 

In  publication  [2]  ,  the  correlation  for  the  directional  drift  of  the 
main  BP  maximum  is  presented  in  relation  to  ♦*  “jjr -f’2  (*,— *oo*ot.  .  For 
the  purposes  of  this  work,  recording  in  relation  to  u  is  preferable 
since  in  the  majority  of  cases  the  length  of  the  antenna  is  considered 
constant. 
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Then  the  decrease  in  dg  (directive  gain) 

.  _  ,  *(0.  a)  _  JLf-.  <0.  m. 

A  -1_"fiioror  2i*  00 1  ' 


It  is  considered  in  (3)  that  (0.  0)=0.  The  ipne  above  designates 

the  operation  of  averaging  for  the  group  of  realizations. 

As  is  apparent  from  (1,2,3),  for  further  computations  it  is 
necessary  to  know  the  statistical  characteristics  of  phase  fluctuations 
in  the  antenna. 


Statistics  of  Phase  Fluctuations  in  the  STWA 

Assume  that  the  technology  for  fabricating  the  original  group  of 
sections  is  such  that  the  wave  number  in  each  section  differs  from  the 
theoretical  k  by  a  random  quantity  Ak.  Then  the  fluctuations  in  the 
phase  in  the  n-th  section  of  the  formed  STWA  ar?  described  by  the 
expression  n 

*.(!/)=/2!**<+0’5fl'c«/^+,)+56Y',  (4) 

i-i  <~i 


where  iK/=*!L  --fluctuations  in  wave  number  in  i-th  section  of  formed 

o 

STWA  classified  as  the  quantity  a 
&Yi  =  ~  --phase  error  at  end  of  i-th  section  governed  by  imper¬ 
fection  in  the  section  junctions,  and  also  classified 
as  the  quantity  o. 

Separate  realization  of  the  phase  fluctuations  in  distribution  is 
shown  in  fig.  la. 

It  follows  from  (4)  that  the  total  errors  #n<y)  increase  with  the 
number  of  the  section,  i.e.  ,  the  quantity  of  the  phase  error  in  the 
n-th  section  is  significantly  affected  by  the  phase  errors  in  the  pre¬ 
vious  sections. 

We  assume  that  the  average  value  Ay i *0 ,  then 
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A  k,  Ay/  =  o,  Ay<Ay /  =  ^  /’ 


(5) 


& 


#•/ 


i  ::  ji 


Figure  1. 

The  statistical  characteristics  Ak^  depend  on  the  statistics  of 
the  wave  number  of  the  original  group  of  sections,  methods  of  rejection 
and  sectioning. 

Assume  that  the  wave  number  of  the  original  group  of  sections  is 

distributed  according  to  the  normal  law  with  zero  average  value  Ak=0 
o 

and  dispersion  oq. 

We  will  consider  the  sections  suitable  for  assembly  of  the  antenna 
if  the  phase  incursion  on  the  length  of  the  section  (equal  to  Akl)  does 
not  exceed  the  tolerance  a,  i.e.,  |Akl|  <_  a.  The  sections  that  do 
not  satisfy  this  condition  are  rejected. 

With  the  selected  method  of  rejection,  the  wave  number  of  the 
"working"  group  of  sections  is  subordinate  to  a  limited  normal  law  of 
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Figure  2. 

distribution: 

-(~5) 

•Si'S#©  1  ’  (6) 

o.  !-*!>  -7-  ♦ 


e"^*dt-- integral  of  probability, 

— parameter  associated  with  the  percentage 
--.y— of  section  rejection. 

r  2o»* 

Thus,  for  example,  1%  rejection  corresponds  to  the  probability 
p(|x|  <_  j)  “  *  (0=0.99,  from  which  £  =1.83. 

The  dependence  of  the  percentage  of  rejection  on  £  is  graphically 
presented  in  fig.  2  (curve  h) .  It  is  apparent  from  the  figure  that 
the  smaller  values  of  rejection  percentage  correspond  to  the  larger  £, 
and  with  £  >  2,  the  percentage  of  rejection  changes  slightly  and  the 
rejection  of  sections  may  not  be  done.  Since  £  is  unequivocally  linked 
to  the  percentage  of  rejection,  the  quantity  £  characterizes  the  degree 
of  economy  of  production. 

We  will  now  examine  the  effect  of  the  sectioning  methods  on  the 

p 

statistic,  of  th<-  ave  number.  For  this  purpose  we  will  examine  the 
selective  (T  .netnod)  and  random  (II  method)  layouts  of  the  antenna  from 


where 


*(!)-=—  l 

V*  .* 


the  "working”  group  of  sections  [2].  With  random  arrangement,  the 
sections  of  the  working  group  are  selected  according  to  the  following 
law.  In  this  case,  the  wave  number  in  the  i-th  section  of  the  STWA 
will  evidently  be  distributed  according  to  the  limited  normal  law  (6). 

With  selective  arrangement ,  sections  are  installed  in  the  first 
place  in  the  antenna  with  positive  phase  incursion  Akl  >  0,  the  second 
place,  with  negative  incursion  Akl  <  0,  etc.  With  this  method  of 
sectioning,  the  wave  number  in  the  i-th  section  of  the  formed  TWA  is 
distributed  according  to  the  limited  one-sided  normal  law. 

By  knowing  the  laws  of  distribution,  it  is  easy  to  obtain  expressions 

for  the  average  value  Ak .  and  dispersion  o  of  the  wave  number  in  the 

i-th  section  of  the  STWA.  These  expressions  are  Dresented  in  [2].  Here 

we  will  only  note  that  for  the  I  method  ,  ZkT  is  constant  for  each 

section,  positive  in  the  odd,  and  negative  in  the  even  sections,  and 

for  the  II  method  Ak^=0 .  Dispersion  o  does  not  depend  on  the  number 

of  the  section  and  is  determined  by  the  dispersion  of  the  original 

2 

group  of  sections  on  and  the  parameter  of  rejection  f; .  The  dependence 
2  2  u 

of  a  /oq  on  this  parameter  for  both  methods  of  assembly  is  presented 
in  fig.  2  (curves  b  and  c  respectively) .  The  expediency  of  using  the 
I  method  of  sectioning  is  apparent  from  the  figure,  since  in  this  case, 
o  is  x'oughly  three  times  smaller  than  in  the  II  method. 


Decrease  in  DG^,  Directional  Dispersion  of  the  Main  Maximum 

Using  different  design  measures  ,  we  will  "  sampld'  the  phase  error 
in  the  antenna  aperture  to  zero  through  M  sections  (fig.  lb).  For 
simplicity  we  assume  that  the  TWA  consists  of  a  whole  number  of  groups 
of  M  sections  each,  i.e.  ,  N=uM.  Then  the  correlations  for  F(0,  o)  and 
Au  look  like: 


F{ 0.  o)=.-l- 

m  M  I 


9  «  I 

SSI 

V—  1  —l 


Iff*  „  u  (») 
e  M  dy 


(7) 


A  u  =  -  JJ  JJ  J  [2  (p v  M - Af)  +  y- N -  1] 0p+(v_„  M 


(y)dy. 


v— I  p—l  —I 


(8) 


Since  errors  at  the  end  of  each  group  are  sampled  to  zero  and  in 
each  group  the  same  algorithms  are  used  for  installing  sections  made 
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of  a  "working"  group  of  large  volume,  then,  evidently,  the  phase  fluc¬ 
tuations  in  each  group  are  statistically  equivalent  and  independent  of 
each  other. 

This  makes  it  possible  to  consider  that 


*0 - gras - M = j  *°p  («/>  'K  v  * *‘- 


where  designates  the  centered  quantities.  The  calculations  of  the 
amount  of  decrease  in  dg  according  to  formula  (3)  with  regard  for  ratio 
(7)  ,  the  average  value  and  dispersion  of  the  course  of  direction  of 
the  main  maximum  using  formula  (8)  are  fairly  cumbersome,  however  do 
not  involve  any  basic  difficulties.  In  the  computations,  it  is  first 
of  all  necessary  to  take  into  consideration  ratio  (9),  then  ratios(A) 
and  (5)  ,  and  also  to  additionally  hypothesize  that  the  correlation 
coefficient  (with  i^j)  of  the  fluctuations  in  the  wave  number  Ak^  equal 
zero.  m  m-i  m  m  o-i 


M 

M- 1 

M 

Af  (?-l 

V  * 

V 

V 

V  V 

P<9 

1 

«--=P+  1 

after  the  appropriate  computations  ,  we  can  write  the  following  final 
results : 

1)  decrease  in  dg  with  sampling  of  the  phase  error  to  zero  through 
M  sections 

4«-sr  )- 

[(GAI  — 4,v — 3)  —  (4.Y  3.1  f)  ^’]  ;  (10) 

2)  average  direction  of  main  maximum 


Au  =  —  aK 


I  A  K,  I  /. 


- M  2AIA’  ’  ' 

a  _  (o  with  M  even, 
where  a^  with  M  odd; 

3)  dispersion  of  the  direction  of  the  main  maximum  of  the  BP 

4.  -  i5;r  [«•  ('»“•- *  (l!r  +  *)~l£  C*’  - 

~o515MMw‘-0(3A<-1)  +  3]}.  (>2) 


Analysis  of  Results 

We  initially  will  analyze  the  effect  on  the  quantities  AJt ,  Am,  o ^ 
of  the  errors  in  wave  number  Ak^  and  phase  errors  at  the  end  of  the 
sections  hy  ^  in  the  absence  of  sampling  of  the  total  phase  error  4>n (y )  - 
We  will  then  evaluate  the  effect  on  the  indicated  parameters  of  a 
sampling  with  =  0  ,  and  finally  ,  we  will  examine  a  general  case  where 
sampling  of  errors  through  M  sections  is  stipulated  in  the  antenna  and 
errors  Ak^  and  Ay^  are  present.  We  will  examine  several  possible  situ¬ 
ations  . 

1.  Sampling  of  the  phase  error  is  not  done  (M=N  ,  u>=l).  In  this 
case,  with  0^=0,  the  ratios  for  AN  ,  Au  and  o^N  have  been  analyzed  in 
publication  [2].  It  is  noted  there,  that  both  for  I  and  for  II  methods 
of  sectioning  ,  with  oL=const ,  distortions  in  the  beam  pattern  are  re¬ 
duced  with  an  increase  in  the  number  of  sections.  This  result  is  almost 
obvious  ,  since  with  an  increase  in  the  number  of  sections  with  L=const , 
dispersion  of  the  phase  errors  diminishes.  When  there  are  errors  which 
arise  at  the  sites  of  junction  of  the  sections,  the  situation  changes: 
increase  in  the  number  of  sections  ,  on  the  one  hand  results  in  a 
decrease  in  the  dispersion  of  phase  errors  governed  by  fluctuations 
in  the  wave  number  ,  and  on  the  other  hand  ,  in  an  increase  in  the  phase 
errors  over  the  length  of  the  antenna  which  are  governed  by  inhomo¬ 
geneities  at  the  sites  of  section  contact.* 

A  similar  situation  must  result  in  the  presence  of  certain  values 
2 

N  at  which  aN  and  ou^  will  be  the  minimum,  i.e.  ,  in  the  existence  of 
an  optimal  number  of  sections  from  which  the  STWA  is  assembled.  With 
fairly  large  N,  where  in  formulas  (10)  and  (12)  with  M=N  and  w=l,  one 
can  ignore  the  terms  of  a  higher  order  of  smallness  than  1 /N ,  the 
minimum  values  for  decrease  in  dg  and  dispersion  of  direction  of  the 
main  BP  maximum  are  reached  practically  simufca  neously  with 

N»  — ,  (13) 

« 


By  using  (A)  ,  it  is  easy  to  show  that  the  maximum  disoersion  of  the 
total  phase  errors  is  reached  at  the  end  of  the  antenna  and  equals: 

«  «*£*  The  first  term  is  associated  with  errors  in  the 

N  •  (ft  '("'Y  . 

wave  number,  while  the  second  is  associated  with  phase  errors  at  the 
sites  of  junction  of  the  sections. 


-Q- 


•  j»»y— r 


where  --ratio  of  standard  value  of  phase  error  at  the  sites 

of  junction  of  the  sections  to  'the  standard  phase  incursion  at  length 
L  in  the  presence  of  fluctuations  in  the  phase  of  the  wave  number. 

Since  ratio  (13)  was  obtained  on  the  hypothesis  that  N  »  1,  then  when 
it  is  used,  the  condition  cL  »  must  be  fulfilled. 

When  this  condition  is  fulfilled  (in  addition  to  the  condition  for 

2 

smallness  of  the  errors  (see  the  footnote  on  p.  00  [sic])  in  which 

the  correlations  (10  and  (12)  were  obtained  .from  formula  (13)  one  can 

determine  the  number  of  sections  into  which  the  antenna  needs  to  be 

"divided”  in  order  to  obtain  the  minimum  distortions  of  the  beam  pattern. 

2 

Since  the  quantity  o  during  selective  layout  is  smaller  than  during 
random  layout  ,  then  the  STWA  can  be  formed  from  a  smaller  number  of 
sections. 

A  o1 

The  dependence  of  the  parameters  rsjr  and  on  the  number  of 

sections  N  with  different  values  of  a  is  presented  in  figs.  3  and  4. 
Figure  3  illustrates  the  dependence  of  the  ratio  ~~  on  the  number 
of  sections  with  random  layout  of  the  antenna  (i.e.  1  Ak,-|  —0) . 

It  is  apparent  from  the  figure  that  in  the  region  of  optimal  N, 
the  quantities  A^  and  change  slowly.  With  a=0,  the  optimums  are 
missing  and  the  distortions  of  the  beam  pattern  are  smaller  ,  the 
greater  the  number  of  sections  N. 

Above  we  examined  a  case  where  the  length  of  an  antenna  made  of 
a  different  number  of  sections  remained  constant ,  while  the  length  of 
the  sections  changed  1=L/N.  When  the  antennais  assembled  from  sections 
of  a  fixed  size  (l=const) ,  the  increase  in  N  results  in  an  increase  in 
errors  governed  both  by  fluctuations  in  the  wave  number  Ak^  and  by 
phase  errors  at  the  sites  of  junction  Ay^.  In  this  case,  the  distortions 
in  the  beam  pattern  increase:  Ajq  increases,  and  the  dispersion  of  the 
direction  of  the  main  BP  maximum  rises. 

2.  A  sampling  is  made  in  every  other  M  section  of  the  phase 

N  2 

errors  (M=  — ) ,  errors  at  the  sites  of  junction  are  missing  (o  =0).  In 

this  case  we  study  the  effect  of  the  "frequency"  of  sampling  the  errors 

on  the  STWA  characteristics.  For  this  purpose,  we  analyze  expressions 

(10),  (11)  and  (12)  with  o*  «0. 
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Figure  3 . 

We  will  first  examine  the  average  value  and  dispersion  of  the 
direction  of  the  main  BP  maximum. 

With  random  layout  of  the  sections  and  (or)  an  even  number  of 
sections  in  the  group  ,  the  direction  of  the  main  BP  maximum  is  maintained 
on  the  average  the  same  as  in  the  absence  of  errors.  When  errors  are 
sampled  after  each  section,  the  quantity  Tu  is  the  maximum.  With  an 
increase  in  the  number  of  sections  in  the  antenna  (with  L=const)  or  in 
the  group  in  the  I  method  of  assembly,  Au  is  diminished.  These  relation¬ 
ships  have  a  simple  physical  interpretation,  if  we  examine  the  phase 
distribution  in  the  group  of  sections.  Thus,  for  example,  with  M=1 
and  the  I  method  of  assembly,  the  distribution  of  phase  errors  in  the 
antenna  is  shown  in  fig.  lc  ,  and  the  average  distribution  of  errors  is 
the  same.  Evidently,  Tu  in  this  case  will  be  the  maximum. 
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Figure  4. 

For  simplicity,  we  will  introduce  the  following  quantity  into  the 
examination 


where  alM  --dispersions  of  direction  in  the  main  BP  maxinun  of 

the  antenna  without  sampling  of  the  phase  errors  and  antenna  with 

sampling  of  the  errors  in  M  sections  respectively.  For  the  examined 

o 

case,  the  quantity  ou^  is  determined  by  the  correlation  (12)  with  N=M 
and  *  0  ,  and  is  analyzed  in  fair  detail,  as  previously  noted,  in 

publication  [2].  Since  then  the  expression  xu  =  n *v.'o^ 

\ 

characterizes  the  gain  in  dispersion  of  the  direction  of  the  main 

maximum  because  of  sampling  the  phase  errors.  As  is  apparent  from  (12), 

N 

^  depends  on  N ,  M  and  the  ratio  between  them  u>=  . 

The  indicated  relationships  are  illustrated  by  the  graphs  in  figs. 

5  and  6  with  a  =  0.  It  is  apparent  from  fig.  5  that  in  the  antenna 

2 

with  constant  number  of  sections  (N=const)  ,  the  greatest' gain  in 
is  obtained  with  sampling  the  phase  errors  at  the  end  of  each  section 
(with  N=12  ,  M=1  ,  Xu~70). 
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With  an  increase  in  the  number  of  sections  M  in  the  group,  the 
quantity  is  diminished  and  the  rate  of  its  decrease  slows  down. 

A  more  significant  "gain'1  in  directional  dispersion  of  the  main  BP 
maximum  is  obtained  in  the  antennas  with  greater  N  (fig.  6,  curves  M= 
const).  However,  with  an  increase  in  M  (reduction  in  u>)  ,  the  dependence 
of  Xu  on  N  becomes  weaker  and  the  gain  practically  becomes  the  same  for 
any  N  with  u=2  (fig.  6).  It  is  also  apparent  from  this  figure,  that 
the  quantity  x  mainly  depends  on  the  number  of  groups  and  its  evalua- 
tion  can  be  made  according  to  the  formula  yulS£~f  ‘ 

In  analyzing  the  decrease  in  dg ,  by  analogy  with  x^,  we  introduce 
into  the  examination  a  quantity  which  equals  the  ratio  of  decrease 
in  dg  in  the  absence  of  sampling  of  phase  errors  to  the  decrease  in 
dg  in  the  STWA  with  sampling  of  errors  through  M  sections  xA=  — 

{see  formula  (10)].  This  quantity  characterizes  the  "gain"  in  reducing 
the  dg  through  the  sampling.  The  dependence  of  on  the  parameters 
of  sectioning  is  illustrated  by  a  series  of  graphs  presented  in  figs. 

7  and  8.  The  numbers  on  the  arrow  in  fig.  7  show  the  quantity  X^  with 
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Graphs  with  a  =  0  correspond  to  the  examined  case 


Figure  6, 


Figure  7 . 


It  is  apparent  from  the  graphs  that  the  nature  of  dependence  of 
06  on  M,  N  and  u  is  mainly  the  same  as  for  Xu.  The  difference  is  that 
with  N  >>  1  ,  the  gain  in  dg  because  of  error  sampling  can  be  estimated 
by  the  formula  x*  «0,3u>+0,4  (in  contrast  to  —)• 

3.  Sampling  is  made  of  the  phase  error  through  M  sections  in  the 
presence  of  fluctuations  in  the  wave  number  and  errors  in  the  sites 
of  junction  (a^  j4  0,  M/N) . 

In  this  case  ,  the  decrease  in  dg  and  dispersion  of  the  course 
of  the  main  maximum  is  described  by  expressions  (10)  and  (12)  respec¬ 
tively. 

As  follows  from  these  expressions  ,  with  M=const  and  with  unchanged 
length  of  the  sections  1  ,with  an  increase  in  the  number  of  sections  N 
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(and  this  means  an  increase  in  the  length  of  the  antenna) ,  the  decrease 
in  dg  remains  practically  unchanged,  and  the  directional  dispersion  of 
the  main  maximum  diminishes. 

Thus,  wwhen  sampling  of  the  phase  errors  is  used,  there  are  no 
limitations  which  are  usually  placed  on  the. length  of  the  antenna  be¬ 
cause  of  the  presenc  of  errors  in  the  APD. 


The  dependences  of  the  quantities  *u  and  on  the  sectioning 
parameters  and  the  size  of  the  errors  at  the  sites  of  junction  a  are 
presented  in  figures  5-8. 

As  is  apparent  from  the  figures  ,  the  appearance  of  the  dependences 
of  gain  in  the  reduction  in  dg  and  dispersion  of  direction  of  the  main 
maximum  with  ajK)  on  the  quantities  u,  N,  M  are  preserved  the  same  as 
in  the  case  of  a=Q.  The  most  significant  discrepancies  in  the  numerical 
values  it  and  x,  are  observed  with  small  M.  This  is  associated  with 
the  fact  that  with  small  M,  the  effect  of  on  the  antenna  parameters 
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are  practically  excluded. 

We  note  that  in  evaluating  the  effect  of  the  sectioning  methods 
and  rejection  of  sections  on  the  decrease  in  dg  and  the  amount  of 
directional  dispersion  of  the  main  BP  maximum,  it  is  necessary  to  use 
the  ratios  (10)- (12)  and  the  graphs  in  fig.  2.  It  is  evident  that  during 
selective  arrangement  of  the  sections  ,  the  quantity  is  smaller  than 
during  random  arrangement.  The  increase  in  the  percentage  of  rejection 
of  the  sections  also  results  in  a  decrease  in  o^j.  Similar  results 
can  be  drawn  in  relation  to  the  quantiy  Aj^. 

In  conclusion,  the  authors  express  their  gratitude  to  Professor 
Ya.  S.  Shifrin  for  attention  to  this  work. 

Appendix 

In  order  to  obtain  a  formula  for  the  course  of  the  direction  of 
the  main  BP  maximum,  we  initially  present  a  technique  to  solve  equations 
of  the  type: 

AuF(u.  fc)  =  0,  fA.n 


where 

F(u,b)--beam  pattern  of  antenna  in  presence  of  errors  in  APD , 

Ascertain  linear  (in  relation  to  u)  operator, 
b--certain  statistical  small  parameter  of  errors  in  APD. 

In  this  case,  one  can  view  (A.l)  as  an  equation  which  implicitly 
assigns  u  as  the  function  of  b,  and  presents  u  in  the  form  of  a  Maclaurin 
series,  limited  only  by  the  first  m  terms: 


*»0 


(A. 2) 


The  derivatives  of  u£(0)  which  are  included  in  (A. 2)  are  defined 
from  the  known  rules  for  differentiating  implicit  functions  [A]. 

Thus ,  for  example  , 


«;  (0) 


Atf  t iu<  W 


») 

oKw+Cfr-  o>K <°>r  i 

A^F H  0) 
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The  quantity  Uq  is  a  solution  to  equation  AuF(u,0)=0,  where  F(u,0) 
is  the  beam  pattern  of  the  antenna  in  the  absence  of  errors  ,  and  is 
usually  considered  to  be  known. 

With  fairly  small  values  in  the  expansion  of  (A. 2),  we  can  limit 
ourselves  to  the  first  two  terms.  Then 

0) 

•  (A.  3) 

With  definite-  ,  f  rom  formula  (A.  3)  one  can  find  the  quantity  of 
expansion  of  the  BP  position  of  zeros  or  equisignal  direction  in  the 
BP,  etc.  In  particular,  with  —  ,  formula  (A.  3)  determines  the 

nu 

directional  course  of  the  main  BP  maximum. 
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STATISTICS  OF  THE  FIELD  OF  FOCUSED  SYSTEMS 


A.  S.  Yevseyev 


Summary 


The  average  characteristics  of  a  focal  spot  are  examined 
in  the  focusing  area  when  there  are  random  phase  errors.  The  cal¬ 
culated  correlations  and  graphs  make  it  possible  to  estimate  the 
decrease  in  field  intensity  in  the  focus  and  the  expansion  of  the 
focal  spot  depending  on  the  numerical  characteristics  of  the  phase 
errors. 


Introduction 

By  focused  antennas  we  usually  mean  those  antennas  in  which  the 
quadratic  component  for  the  differene  in  the  course  of  beams  from  the 
antenna  to  the  observation  point  located  in  Fresnel's  zone  is  com¬ 
pensated  for  in  some  manner. 

It  is  common  knowledge  that  in  Fresnel's  zone,  the  expression  for 
the  difference  in  the  course  of  the  beams  (fig.  1) 

«(#-*,)*« sinO-  +  ...  (i) 

does  not  take  into  consideration  the  terms  of  the  series  above  the 
second  degree  in  relation  to  z.  By  artificially  introducing  phase 
distribution  of  the  field  along  the  antenna  equal  to  .  we 

can,  in  the  limits  of  small  angles  in  the  region  9 m  8q  at  distance  Rq 
(focal  distance) ,  obtain  the  same  field  as  in  the  distant  zone  of  the 
antenna.  With  R^Rq  ,  the  distribution  of  the  field  is  different  than 
in  the  distant  zone. 


Figure  1. 

Interest  in  the  focused  systems  has  recently  noticeably  increased 
since  they  can  be  used  to  solve  a  number  of  important  problems.  Thus, 
for  example  ,  the  use  of  focusing  permits  measurement  of  the  beam 
pattern  of  large  antennas  at  reduced  distances  ,  expansion  of  the 
potentialities  of  the  antennas  of  side- scanning  radar  stations,  etc. 

In  this  respect,  it  is  definitely  important  to  study  the  para¬ 
meters  of  the  focused  systems  when  there  are  inevitably  existing 
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different  destabilizing  factors.  This  work  examines  the  average 
intensity  and  width  of  the  focal  spot  of  a  focused  antenna  when  there 
are  phase  errors  along  the  antenna.  Phase  errors  can  be  dictated  both 
by  inaccurate  fabrication  and  instabilities  in  the  parameters  of  the 
elements  of  the  antenna  itself  (internal  errors)  ,  and  by  the  presence 
of  fluctuations  in  the  refraction  index  of  the  medium  (external  errors). 


Original  Correlations 

The  field  intensity  of  the  focused  linear  antenna  in  Fresnel's 
zone  with  accuracy  to  a  constant  multiplier  equals  [1]: 

P*W<  a)  =  |j^(.v)/l(.Vi)eiJ,A  x,Hla(*i ~^dxdxlt  (2) 


2z 

where  A(x) --amplitude  distribution,  x=  £—  --relative  coordinate, 

^  (sinfl— sin00)  --generalized  angle,  ^--direction  of  the  radiation 

.  _  kL*  /cos5  0  cos4  tl„  c 

maximum  ,  a=— •(  — - - — |—  focusing  parameter. 

o  \  R  R  o 


In  the  case  of  uniform  amplitude  distribution  A(x)=Aq=%,  the 
expression  for  PQ(ifi  ,a)  looks  like 

P.  (*,  a)  =  {[C  (t)  +  C  (S)P  +  (S  (t)  -f  S  (|)JS),  (3) 

where  t -a(l— *  =  0(1  +  —  ';*.  C(t).  S(t)—  --Fresnel's  integrals. 

With  R=Rq  and  g^Sq,  one  can  consider  that  a=0.  Then  P„ (tf,  0)  -  . 

The  transverse  dimensions  of  the  focal  spot  (at  the  level  of  half 
power)  correspond  to  the  width  of  the  antenna  BP  in  Fraunhofer's  zone, 
and  are  determined  by  the  known  correlation:  2 -,,t'  With  i(/=0, 

we  have  t  = £  =  a  ,  and 

P .  W>,  a)  =  P,  (0,  a)  —  —  [C  (a)  4-  S*  (a)]. 

a 


The  graph  for  the  dependence  of  Pq(0,3)  on  the  focusing  parameter  is 
presented  in  fig.  2  (solid  curve).  By  assuming  PQ(0,a)=0.5,  we  find 
the  longitudinal  dimensions  of  the  focal  spot 

The  ratio  of  the  longitudinal  spot  dimension  to  its  widtn  in  relative 

units  x i ,(*. 

Mo.sr  ~ 
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Average  Field  Intensity 

When  there  are  random  phase  errors  <j>(x)  along  the  antenna,  the 
individual  realization  of  the  field  intensity  will  be 


=  _L  jj e<  C*  <*>-»  <*.)!  e‘  ♦  <*-*'»+'»  (  ■*') 


P(yp,  a) 


dxdxx. 


The  statistics  P(i|»,a)  depends  on  the  numerical  characteristics  of  the 
phase  errors  <j>(x).  We  will  examine  two  types  of  errors:  internal 
and  external. 

The  following  assumptions  are  usually  correct  for  the  first  type 
of  errors:  the  errors  are  small  and  stationary  over  the  longth  of  the 
antenna,  i.e.  ,  the  average  value  if  (x)  =0  ,  the  dispersion 


and  the  correlation  coefficient  r(x ,x1> =r (x-x1> .  It  evidently  makes 
sense  to  study  the  effect  of  these  errors  on  the  field  intensity  for 
such  focusing  areas  where  the  standard  value  of  the  phase  considerably 
exceeds  the  quantity  of  the  negligible  third  term  in  the  series  in 


expression  (1),  i.e., 


kL%  »ln  9 
16  R* 


At  the  same  time  [1] 


*>0,621  L*ll. 


The  latter  inequality  determines  the  near  boundary  of  the  Fresnel's 
zone.  When  the  wave  spreads  in  the  turbulent  medium  ,  the  errors 
(fluctuations)  may  have  a  random  quantity.  To  describe  the  errors  , 
the  literature  uses  both  Gaussian  and  Kolmogoroff  models  of  an  inhomo¬ 
geneous  medium.  We  will  examine  both  models  below. 
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1.  We  will  first  study  the  average  field  intensity  when  there 
are  errors  of  internal  origin.  .By  averaging  the  correlation  (4)  for 
a  group  of  realizations,  or  for  time,  we  obtain 

P <*/  a)  =  P.  W.fa)  ( 1  - a*)  4-  o'  Y  (>f.  a). 


where 


+i 


Y  W>,  a)  -  ~  1 1  r  (*—*,)  e  dxdxx. 


(5) 


-i 


We  will  further  adopt  the  correlation  coefficient  of  phase 
errors  in  the  Gaussian  form  r=exp|  —  t  ,  where  c=  ~~ 

relative  radius  of  the  error  correlation.  The  integral  Y(i|/  ,a)  for 

this  case  is  computed  aoproximately  in  appendix  1  (the  computation 

used  the  results  of  publication  [33). 

Formula  (5)  and  the  results  of  appendix  1  make  it  possible  to 

calculate  the  average  field  intensity  with  different  values  of  the 

2 

parameters  a ,  ,  c  and  a  .  , 

With  a  =iO,  the  expression  for  P(^  ,0)  coincides  with  the 
analogous  expression  for  the  average  beam  pattern  in  the  distant  zone 
that  was  studied  in  detail  in  publication  [2].  It  also  presents  curves 
which  characterize  the  decrease  in  field  intensity  in  the  main  direction, 
and  curves  which  determine  the  width  of  the  beam  pattern.  These  curves 
make  it  possible  to  determine  the  decrease  in  field  intensity  in  the 
focus  (fig.  2.3)  of  work  [2])and  the  transverse  dimensions  of  the  focal 
spot  (fig.  2.15,  work  [2]). 

With  *=0, 


P (0,  o)  =  P,  (0,  a)  (1  -a*)  +  o*  Y  (0,  a),  ( 6) 


where 


V(0,  a) 


■tJJ 

-1 


dxdxv 


Figure  2a  ,b  shows  the  dependence  of  F(0 ,a)  on  a  for  a  series  of  values 
o 

of  o  and  c.  The  curves  of  fig.  3  define  the  width  of  the  focal  spot 
in  a  longitudinal  direction. 
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Figure  3 . 


Figure  A . 


By  analyzing  figs.  2,3  and  A  we  can  draw  the  following  conclu¬ 
sions  . 

a)  The  presence  of  phase  errors  results  in  a  decrease  in  field 
intensity  in  the  focus  and  expansion  of  the  focal  spot  ,  ana  the  trans¬ 
verse  and  longitudinal  dimensions  of  the  spot  ip/"- ease  roughly  equally 

0  2Afl05p  ’  \  ' 

(thus,  for  example,  with  c=0.5  and  o  =0.5  — — 1 — *=  i,65  ; 

z&S’o.sp  1 

b)  Reduction  in  the  relative  radius  of  the  correlation  of  phase 
errors  or  increase  in  their  dispersion  causes  expansion  of  the  focal 
spot  and  reduction  in  field  intensity  near  the  focusing  point. 

2.  We  will  now  examine  the  effect  of  turbulence  of  the  medium 


on  field  intensity  of  the  focused  systems  ,  after  limiting  ourselves 
to  consideration  of  only  phase  fluctuations.  If  in  order  to  describe 
the  medium,  we  adopt  the  Gaussian  model  of  an  inhomogeneous  medium, 
then  in  this  case,  the  results  of  the  calculations  obtained  in  section 
1  will  be  completely  correct.  In  this  case  one  should  only  bear  in 


mind  that  the  dispersion  of  phase  errors  is  determined  by  the  ratio  [2] 

oj_0!1  {\n)*K*kR,  where  T5n> ^--dispersion  of  fluctuations  in  the 

2 

refraction  index  of  the  medium,  while  the  radius  of  the  correlation  c 


is  roughly  equal  to  the  average  scale  of  inhomogeneity  1^. 

The  calculation  looks  somewhat  different  for  the  case  where  the 


Kolmogoroff  model  of  turbulent  medium  is  used. 
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In  this  case  (4] 


t 

e'  ('  <«>-•  wi  =  exp  l  -Cl  I  Jf ~ X,  |S/J] . 
where  4  ' 

Here  R-- length  of  wave  path  in  turbulent  medium;  B- -parameter 
which  characterizes  the  intensity  of  the  dielectric  inhomogeneities 
(B  ^(lO'^-lO"^)  cm-^^]  ,  on  the  average  B=1.3  x  10-1^  cm’^^ 

j  «ithi*-*il>»/r,-/?- 

The  5/3  lav;  for  the  structural  function  of  the  phase  is  confirmed 
all  the  way  to  the  bases  jz-z^l^  km  [5]. 

As  previously,  in  the  absence  of  fluctuations  in  the  phase,  the 
field  intensity  is  determined  by  the  ratio  (3). 


Figure  5.  Figure  6. 


The  average  field  intensity 


/  />(*,  a)  =  -M)e  ’  '  1  'dxdx v 


The  integral  in  correlation  (7)  is  computed  approximately  in  appendix  2. 
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Figure  8. 


By  assuming  a=0 ,  we  find  from  correlation  (7) 

_  .  jfi  —C*  I*— (*—*,) 

mill-* 


/>».  0)_- 

—I 


dxdxv 


The  curves  for  average  field  intensity  which  are  computed  from 

2 

this  formula  for  the  two  values  are  presented  in  fig.  5.  The  graph 

(fig.  6)  characterizes  the  dependence  of  field  intensity  near  the 

focusing  point  on  the  parameter  a.  Figure  7,8  presents  the  graphs  which 

2 

determine  the  effect  of  on  the  longitudinal  and  tranverse  dimen¬ 
sions  of  the  focal  spot.  These  graphs  can  be  used  to  estimate  the 
linear  dimensions  of  the  spot  at  the  assigned  distance  when  the  antenna 
is  transmitting.  From  an  analysis  of  the  presented  graphs  we  can  draw 
the  conclusion  that  with  an  increase  in  the  phase  errors ,  the  field 
intensity  drops  ,  while  the  transverse  and  longitudinal  dimensions  of 
the  focal  spot  increase  equally. 

When  the  graphs  are  used,  it  is  necessary  to  bear  in  mind  that 

2  2 
the  quantities  and  a  are  mutually  dependent ,  since  both  and  a 

are  determined  by  the  quantities  R,  L  and  x.  The  distance  Rq  and  the 

patameter  of  turbulence  B  are  independent  parameters. 

In  conclusion,  the  author  expresses  deep  gratitude  to  Ya.  S. 

Shifrin  for  supervising  this  work,  and  L.  G.  Korniyenko  for  advice 

and  attention  to  it. 
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Appendix  1_. 

1.  l.Te  cannot  precisely  compute  the  integral 


~  ~t/l>  +H  <*-*,)+!•  (  1|— jr*} 


dxdx t 


with  random  c  and  a.  Therefore  we  will  examine  the  asymptotic 
expressions  for  Y  with  c<<1  and  c>>1: 
a)  with  c<<1 

By  using  the  variables  u  =  and  >  we  ^ave 


0 

I  - 


I'2  >T-«  ,  0  _  yr^a 

\  c  '•  dv  j  e'»u’du+  J  e  <*  '  % 

—  <  *  2" — a)  -  »T 


»T  _  *!_  _]  yftu 
e dv~2  j  e  <* 

—  (  lT-f  a)  0 


x  J 


au 


_/2u 


Passing  further  to  a  new  variable  y~-  -  ,  we  obtain 

Y=2 


tin  Tacy  ( I  —  ~) 


(eu  , 

I  —  — 

- 2—dy. 

J 

»/« 

-4/c4 


With  accuracy  to  the  quantity  — g  e”  '  we  have 


—2: 

S  »'p, 


Vp,  i  Yp,  t  Vp,  7VT, 

j*  F(x\)dr\—  J  F(i])dr\  +  J  F(t|)dtj—  j  F(tj)dt] 


- —  arctg  p, 

o 


(A.  1) 


(A. 2) 


(A.  3) 


where  *l'-*“2ac ±c^  P  =  nc*  *U»)  =*~  \  c  '  dt 

probability. 


--integral  of 
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*rr 


>W.*. 


By  designating  2 *7i  ~*n  0,  +  ,^’,,  2,  rt  ~2**  a*  l  'h*  and  by 

taking  into  account  Pi  =  P\.  &(**)  ,  with  the  help  of  the  known 

expression  for 

F(x)  -y 

'  Zj  13..  .(2n  +  l)<2*  +  2) 


we  will  write  the  integral  in  the  following  form: 


Re  Y  =  I  n(Re  0  (?„)  +  Re  <f>(;„)|  4  JJ  X 

l  H---0 


(-1)"  T  Im#*' 


■  13  .  .  .(2«H- 


— * n- —  f  ->  y. 

l)(2n+2)  U 


(  —  1)"  2"  Ini  7-r' ' 
7T~.  .(2n  -f  l)(2>H  2) 


— 2arctgp!  . 


(A.  4) 


Expression  (A.  3)  has  a  broad  area  of  application,  practically  to  1, 

With  c<_0.5,  by  assuming  &<1  ,  we  have 


ReY*i  {*  [0  j + 0  (f )] f  m  F(f)-hr"pF(f)-2p} 


(A. 5) 


b)  c>>1. 

_x2  2 

By  using  the  expansion  of  the  exponential  function  e  ~  1-x  , 
we  will  write  the  correlation  (A.l)  with  accuracy  to  terms  of  the 
second  order  of  smallness  in  the  form: 


r  i 

Re  V  =  4P, (if,  fl)— —■  J  u cos  V2^usin2au(\  T—u)du. 


(A. 6) 


Here  Pt($,  a) a)F*(\f,  a), 


M«.  a) ix~t~  /i{c[«(l-J)']  + 

+ci«(,+i),j-,si,(i-in-,sK,+»)‘B- 

P,  w.  a)=TI  <lC  w  +  C +  Is <t)  +  S {m‘ 

**  »■ 

The  integral  in  expression  (A. 6)  is  computed  in  appendix  2. 


ix  2 


We  will  examine  the  integral 


-1 


-cj  |*-x,|8/3+l*  (»— Jr,H  M  (  Jtf—Jt*) 


dxdxx. 


(A. 2.1) 


By  replacing  the  variables  t  v==i±Ji  ,  we  transform  (l)into 

2  2 

the  form: 


o 


^~mmdu  j"  e,4e"(to- 


5/3 


-  <1— W> 

— ifrt*  sin4au(l  — n)  ^ 
eu 


i 


(A. 2. 2) 


It  is  not  possible  to  compute  (A. 2. 2)  in  the  quadratures.  For  an 

approximate  computation,  we  will  use  the  method  of  approximating 

functions.  The  most  suitable  is  the  function  r(u)=cos  1.2  u-0.08  x 

o 

sin  n  u  (fig.  9.).  We  will  write  the  quantity  3.18  +  «»  where 

n--the  nearest  whole  number.  Then  I  a  |  <_0 . 5  and  ra»  1-0. 72otu^-aO.  08  sin 
n  u.  It  is  easy  to  see  that  in  the  expression  for  rn+a  there  are 
trigonometric  functions  of  degrees  l,2,...h.  By  using  the  expansion 
of  these  functions  through  the  function  of  multiple  arcs  ,  one  can 
reduce  the  expression  (A. 2. 2)  to  the  sum  of  integrals  of  the  type: 


u)  cos  v  a  , 

—  <lu 

sin  y  u 

After  the  computations  ,  we  have 

i cos(4aB*-a&?)- ~  +  ^-}  (lafl;)j cos 

+  'T  V slno^ilC^)— C(4aB*)|  +  ^cos(-l<i^— aftj)— 
“T  /•^cos^[S(a^)  +  S(4^)l+  ^  1  ^sin,^| C(«iA=)  - 

+  C(4flB‘)j|  , 

lt  ~  4*  {4i sin  (4a®i — a*D  y  l^C0Sfl,,.’lc(^Br)-cK)l  - 

— T  /-£sinafcilS  (4a/>-)— S  sin(4oK-<i^)~ 

“7-  )  5  c°sq6s  IC  (4aBD  +  C (a6Dl~  Y  sin  a<>’  x 

X|S(4aB*)+S(o^)]j, 


l 

/,  j  -- J  u  sin  4<jm  ( I 


■ 

■5 


sin4au(l  — 


where 


6,  „£-l.  fc~i+l,  B.-1  +  i-, 

Mf.-) 
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EFFICIENCY  AND  INPUT  IMPEDANCE  OF  PLANAR 
EMITTER  PLACED  NEAR  THE  INTERFACE 


V.  P.  Kul’tsep 


Summary 


Strict  general  formulas  are  derived  for  complex  input  impedance 
of  a  random  plane-parallel  interface  of  an  emitter  with  assigned 
current  distribution.  A  calculation  is  made  of  the  dependence  of  the 
input  impedance  and  efficiency  of  the  inphase  planar  disk  emitter  on 
its  radius  and  the  distance  from  the  interface  for  emitters  of  any 


Introduction 

When  an  emitter  is  placed  near  the  interface  of  two  media  1  and  2 
with  different  permeabilities  and  C2  .  the  emitted  energy  is  generally 
redistributed  towards  the  medium  with  greater  permeability.  If  the 
radiation  towards  medium  1  is  useful ,  then  the  efficiency  of  the 
emitter  is  defined  as  the  ratio  of  the  active  power  emitted  into  medium 
1  to  the  complete  active  power  fed  to  the  emitter. 

It  is  important  to  determine  the  maximum  efficiency  that  can  be 
obtained  from  a  planar  emitter  of  assigned  finite  parameters  which 
is  placed  parallel  to  the  interface. 

Initial  Correlations 

We  will  use  the  Hansen-Bekerly  method  to  compute  the  efficiency. 

It  has  been  used  repeatedly  to  determine  the  efficiency  of  emitters 
located  near  the  interface  in  the  same  medium  1  to  which  it  is  desirable 
to  direct  the  radiation  maximum  [1,2]. 

We  will  examine  another  case  here  where  the  emitter  is  located 
in  medium  2,  while  it  is  desirable  to  direct  the  radiation  to  medium  1. 

Assume  that  the  emitter  is  assigned  in  the  form  of  distribution 
of  the  current  density  vector  i  in  a  certain  volume  V.  One  can  show 
(see  the  appendix)  that  the  complex  power  which  is  transferred  through 
the  plane  z=z^  is  located  below  the  interface  ,  but  above  the  emitter 
(fig.  1), 


P,  (z)  —  60.i*k  £ 

«»-•  Lo 

v  >*(|  ra.c-^Jc-'"" .  ,  Ifc’l*  ~1  • 

nij  *2  •'  *  * 


The  complex  power  which  is  transferred  through  the  plane  z=z^  lo 
cated  below  the  emitter: 


Pt  («•)  =  60n*  k  J  [  f (  /„  +  «.&)  (  5?-+ 

+ ■ ~  f  ( fc + «.  id  ( fc + . 

*2  0  *  J 
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Figure  1. 


In  these  expressions: 

k,  =.  K  |  k  2n/X,  h-  wavelength 

mt—  |,'r  cj — /*  •  e|*«i+i60ffA-  (3) 

l*i  mt  —  Mtmi  e^m*— e,m, 

a,  - - - ; - .  O* - -  • 

Mi  mi  +  Mi  mi  *i  mi  +  **  *i 

Mi .e.^  ,oi--magnetic  ,  dielectric  permeabilities  and  conductance  of 
media  1  and  2. 


V 


--coefficients  of  expansion  of  the  vector  potential  according  to  the 
system  of  orthogonal  vector  functions  whose  projections  onto  the 
axes  of  the  cyclindrical  coordinate  system 

^-0-  -'fc-f  '.<»#■  <s> 

**1 

where  <f.  z  --coordinates  of  the  points  of  volume  V  in 

which  the  distribution  of  the  current  density  vector  i  is  assigned, 

J  --Besselian  function;  the  asterisk  designates  the  complex-conjugated 
n 

quantity. 

Expressions  (1)  and  (2)  generally  do  not  make  it  possible  to  find 
the  complete  power  fed  to  the  emitter  which  is  located  in  the  medium 
with  losses,  since  they  do  not  take  into  consideration  the  power 
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which  is  stored  and  can  be  scattered  in  the  volume  between  planes 
z=Zj  and  t.=z^  If  however,  we  examine  a  particular  case  of  a  planar 
emitter  with  zero  extension  on  the  z-axis  ,  and  we  assume  that  z^=Z2  = 
-h ,  then  expressions  (1)  and  (2)  make  it  possible  to  find  quite 
accurate  values  for  the  assigned  distribution  of  currents:  Pj(-h)-- 
complete  power  emitted  to  the  upper  half-space;  P2 (-h) --complete 
power  emitted  to  the  lower  half-space  and  P  =P,(-h)  +  P0(-h)-- 

ftX  jL  z 

complete  power  fed  to  the  emitter. 

After  dividing  the  complete  power  fed  to  the  emitter  into  the 
square  of  the  effective  current  value  at  the  emitter  inlet  Iq  ,  we 
find  the  complete  input  impedance  of  the  emitter 

Z*x  -  R,x  f  i  (6) 

If  we  substitute  the  coordinate  of  the  interface  Z=0  into  (1)  , 
then  we  obtain  the  power  P-^(O)  which  passes  through  the  interface  into 
medium  1  ,  and  the  corresponding  impedance 

Zx-«s  +  IXr-P,(0)'/g.  (7) 


By  definition  given  above,  the  efficiency 


We  will  introduce  the  designations 


C„=f  e*'"*'1  Cji  =  ±JL  « 

'•  /,  m,  '*• 


Since  in  our  case,  i2=0,  then 

7-y-(/r)-/^y-,(/0]e-,"V</rrf9, 
,_l_  JJ  JL  Jm(tr)  +  il,XUn_l(lr)]'-,'*rdrd'f.  (10) 


where  iz ,  i^ ,  ir--proj ections  of  vector  i  onto  coordinate  axes,  S-- 
area  occupied  by  the  assigned  current  distribution. 
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In  addition,  we  designate: 

*» 


U'~T  S  M-  u>  -jr  J  M- 


(ii) 


Then  expressions  (6)  and  (7)  for  impedance  of  radiation  and  the 
input  impedance  can  be  written  in  a  form  convenient  for  calculations: 


z, = — —  f  r  — j-  ( ,  —  i  aa  [*)  m*  +  <  i — i  «a  n  c--j— '  . 

*  jJ  |  "«i  hY  I 

!„  -  l^;iO-a,c . •)V,1 ■ 

K  i  I  "h  *i‘  I 


(12» 


WS) 


The  cofactors  of  the  integrands  in  (12)  and  (13)  have  a  specific 

physical  meaning.  Thus,  u  --beam  patterns  of  the  emitter  which  have 

s 

been  averaged  for  the  azimuth  angle  in  a  free  space  for  peroendicularly 
(u2)  and  parallel  (u^)  polarized  waves.  The  multipliers  before  ug 
show  how  the  beam  pattern  of  the  emitter  changes  because  of  the 
presence  of  an  interface.  The  cofactor  ldl/n^  is  the  element  of  solid 
angle  in  medium  2  multiplied  by  k. 

Input  Impedance  of  Disk  Emitter 

Figure  2. 
a 


One  can  assume  that  among  the  planar  inphase  emitters  of 
assigned  area ,  the  emitter  of  round  shape  with  continuous  and  uniform 
distribution  of  the  sources  has  an  efficiency  close  to  the  highest. 

We  will  examine  the  planar  current  disk  of  radius  L  (fig.  2) 
which  is  located  in  medium  2  parallel  to  the  interface  at  depth  h. 

It  has  a  distribution  of  the  current  density  vector  which  is  constant 
in  quantity  and  direction.  Its  projections  look  like: 

<? 

/,  =/cos^p,  =  / sin q>.  (14) 
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The  direction  of  the  current  density  vector  is  shown  in  fig.  2  by 
the  arrows. 

Current  density  I  is  linked  to  the  input  current  Iq  by  a  coeffi¬ 
cient  which  depends  on  the  plan  for  connecting  the  source  to  the 
emitter.  We  will  adopt  this  coefficient  such  that  the  current  density 
on  the  disk  equals  the  current  density  in  a  square  emitter  which  is 
equidimensional  in  area.  Then  I=lQ/L/~ii  .  This  essentially  means 
that  the  power  source  is  connected  to  the  disk  by  chord  aa  whose 
length  equals  the  length  of  a  side  of  the  equidimensional  square  shown 
in  fig.  2  by  the  dotted  line. 

By  substituting  conditions  (14)  into  (10)  in  this  case,  and  by 
integrating  in  the  limits  of  the  disk  area  ,  we  find 


ikJ'm 


for  n=  ±1 


(15) 


for  the  remaining  n. 


By  subst  ituting  (15)  into  (11)  ,  and  further  into  (13)  ,  we  obtain 
an  accurate  expression  for  the  input  impedance  of  the  disk  with  uni¬ 
form  distribution  of  the  current  density  located  near  the  interface 
of  the  medium  with  any  parameters. 

We  will  examine  certain  particular  cases.  Assume  that  the  emitter 
is  located  on  the  interface  (h=0)  of  two  media  without  losses.  Then 
with  1  <  k2  ,  the  integrand  in  (13)  will  be  real,  while  with  1  >  k2  it 
is  imaginary.  By  correspondingly  dividing  the  integration  interval 
in  (13)  into  two  parts,  and  with  1  <  k2 ,  making  a  replacement  with  the 
variable  l=k2sfn  # ,  we  have 


j  tin  j  a 

f  ys ((ClL sin  8)l(l  +  o«)  + cos* » (| -«,)!—• - 

V'*«  \  J0 


where ^2- -angle  between  the  z-axis  and  the  direction  to  the  observation 
point. 

It  is  apparent  from  here  that  the  integrand  in  the  first  integral 
is  the  disk  beam  pattern.  Thus,  the  disk  example  shows  the  known  situ¬ 
ation  where  the  active  part  of  the  input  impedance  of  the  emitter  is 
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determined  by  integrating  its  bean  pattern  in  the  region  of  real  angles. 
It  is  new  that  for  the  selected  power  circuit  .both  the  active  and 
reactive  part  of  the  input  impedance  are  strictly  defined  in  Ohms. 

We  will  find  approximate  expressions  for  the  input  impedance ,  by 
using  the  fact  that  a  large  part  of  the  emitted  power  is  concentrated 
in  the  region  of  the  beam  pattern  maximum,  i.e.  ,  near  1^=0. 

If  the  dimensions  of  the  emitter  are  great  O^L  >>  1)  ,  then  in  the 
first  term  (16)  one  can  assume  that  m2=k2,  m^-k^,  sinl?  =  and  cos1?  = 

1  ,  while  the  upper  limit  of  integration  is  equal  to  infinity.  Then 


I20n 

VT.+  lT, 


(17) 


If  the  emitter  is  located  in  a  uniform  medium  (02=03=0) ,  then  its 
input  impedance  equals  60it//T^  --half  of  the  wave  impedance  of  the 
medium.  On  the  contrary ,  if  the  interface  is  completely  reflecting 
(02=  -03=1)  ,  then  the  input  impedance  of  the  emitter  equals  120n//c"^  . 

As  follows  from  (17)  ,  the  active  part  of  the  input  impedance  of 
a  large  emitter  located  on  the  interface  lies  between  these  extreme 
values. 

If,  on  the  contrary,  the  radius  of  the  emitter  is  very  small, 
then  the  Besselian  function  in  the  first  integral  (16)  can  be  replaced 
by  the  first  term  in  the  series,  and  by  assuming  m2=k2  and  m^=k^,  we 
find 


R 


MX 


•lOn  (k*L)' 


_ l _ 

>7,+  )  T,  ' 


(18) 


Figure  3  shows  the  dependences  of  R  on  the  disk  radius  L  which 
are  obtained  by  numerical  integration  (16)  with  different  values  of  £2* 
It  is  apparent  from  the  curves  that  if  the  radius  of  the  disk  is 
small ,  then  the  input  impedance  rises  in  proportion  to  the  square  of 
the  radius  of  the  disk  in  accordance  with  (18).  With  k2L=1.5,  the 
rise  in  impedance  is  slowed  down,  and  with  k2L=10  ,  the  impedance  prac¬ 
tically  reaches  its  limiting  value,  determined  by  (17). 

The  reactive  part  of  the  input  impedance  is  determined  by  the 
second  integral  in  (16).  It  is  easy  to  see  that  with  1  -*■  *» ,  the  average 
value  of  the  integrand  diminishes  as  Ml.  This  does  not  guarantee 
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-- 


convergence  of  the  integral.  This  physically  means  that  in  order  to 
maintain  rectangular  distribution  of  the  current  ,  it  is  necessary  to 
create  around  the  emitter  an  infinitely  large  reactive  field,  and  con¬ 
sequently,  the  assigned  distribution  many  not  be  realized  with  any 
value  of  L.  If  we  examine  the  medium  with  losses,  then  with  rectangular 
distribution  of  the  current,  the  first  integral  in  (16) differs.  It 
determines  the  active  part  of  the  input  impedance. 

We  will  therefore  examine  current  distribution  which  diminishes 
to  zero  on  the  edges  of  the  emitter,  for  example: 


(19) 


The  current  density  in  the  center  of  the  emitter  is  assigned  here 
to  be  the  same  as  in  the  emitter  examined  above  with  uniform  distribu¬ 
tion  of  the  current  density. 

The  expression  for  the  input  impedance  with  current  distribution 
diminishing  to  the  edge  will  be  similar  to  (16).  It  yields  the  final 
result  both  for  the  active  and  for  the  reactive  parts  of  the  input 
impedance  with  any  parameters  of  the  medium  2. 

As  an  example  we  will  find  one  of  the  components  of  the  active 
part  of  the  input  impedance  associated  with  waves  of  parallel  polariza¬ 
tion.  The  real  part  of  the  second  term  in  (13)  corresponds  to  it. 

The  precise  expression  for  this  component  looks  like: 


J  JUIL)  Rc  (I  -a,)  j  £  . 


(20)  • 
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In  this  case  ,  the  convergence  of  the  integral  is  guaranteed  by 
the  fact  that  the  integrand  with  /  -*•  «  diminishes  as  VP. 

We  will  approximately  compute  part  of  the  integral  in  the  limits 
from  1/L  to  ».  If  the  emitter  is  so  small,  that  |k2L|<<  1,  then  even 
on  the  lower  limit  of  integration,  one  can  assume  that  m2=m-^=il.  In 
addition,  we  replace  the  square  of  the  Besselian  function  by  its 
average  value  for  the  period  1/nlL. 

As  a  result ,  we  have  for  the  semiconducting  medium 


Thus,  the  constant  component  of  the  emitter  losses  is  obtained. 
The  remaining  components  of  the  input  impedance  can  also  be  obtained 
from  the  general  expression  (13)  with  the  help  of  the  approximations 
suggested  in  [2],  or  numerically. 

Input  Impedance  of  Square  Emitter 

We  will  now  examine  the  case  of  a  planar  emitter  of  square  shape 
(fig.  A).  We  will  assume  that  its  area  is  equal  to  the  area  of  the 
round  emitter  with  radius  L  examined  above.  Then  the  length  of  a  side 
of  the  square  emitter  2a=/irL.  We  will  assign  uniform  distribution  of 
the  current  density  of  the  type  (1A). 

Determination  of  the  coefficients  c  with  the  help  of  (10)  in 

Oil 

the  case  of  a  square  emitter  runs  into  certain  difficulties.  We  will 
therefore  use  the  following  concept  for  beam  patterns  of  the  emitter 
which  have  been  averaged  for  the  azimuth: 


ut  =  -~  f 

&»•/’  J 

9  — n 


to,'o  in 


J 


JjV,  cos 5  +  iy  sin  I) eM  ‘""•I-**0**  dxdy 

9 

ff<»  sin  l  +  iy  cos  \)  'dxdy 


dl 


dt 


(22) 


where  A —azimuth  coordinate  of  the  observation  point,  x,  y  .--rectangu¬ 
lar  coordinates  of  the  points  on  the  emitter. 

Expressions  (22)  can  be  obtained  from  (11)  with  the  help  of  a 
generalized  equation  of  closure,  presenting  the  sum  of  the  squares  of 
the  moduli  cgn  in  the  form  of  an  integral  from  the  product  of  the  sum 
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Figure  A. 

er 

V  , 

%  times  the  sum  complexly  conjugate  to  it  ,  and  by  further 

transforming  the  obtained  expressions  with  the  help  of  (10). 

By  assigning  ix=lQ/2a,  iy=0  and  by  integrating  in  (22)  for  the 
area  of  a  square  with  side  2a,  we  find  the  beam  pattern  of  the  square 
emitter  which  is  averaged  for  the  azimuth. 

n/2 

—S55  j  (23) 

-S/2 


If  the  dimensions  of  the  emitter  are  large ,  then  approximately 


it*  sin*  a/ 

2nji«  /*  ’ 


The  expression  for  u^  is  obtained  analogously.  Assume  that  the 

emitter  is  located  on  the  boundary  of  two  media  without  losses. 

Assuming,  as  with  the  derivation  of  (17),  m2=^i  »  an<*  the 

upper  limit  of  integration  is  equal  to  infinity  ,  we  find  that  the 

input  impedance  of  the  square  emitter  is  expressed  by  formula  (17). 

Thus  ,  the  input  impedances  of  the  square  and  round  emitters  of  large 

size  are  equal,  if  their  areas  are  equal.  If  the  square  emitter  is 

2  2  2 

6mall ,  then  U2=u-j=k  a  /2-n  . 

Consequently,  all  the  conclusions  drawn  above  in  relation  to  the 
input  impedance  of  the  round  emitter  are  also  correct  for  the  square 
emitter  of  the  same  area.  We  now  hypothesize  that  the  square 
emitter  is  divided  into  a  certain  number  of  rectangular  elements  with 
correlation  of  the  sides  c/d,  arranged,  as  shown  in  fig.  A,  and  con¬ 
nected  parallel-in  series.  It  is  evident  that  the  input  impedance  of 
the  element  is  c/d-times  larger  than  the  input  impedance  of  the  emitter. 
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If  the  element  has  a  square  shape  ,  then  its  input  impedance  is  the 
linear  impedance  of  the  emitter.  Consequently,  the  graph  in  fig.  3 
yields  not  only  the  input  impedance  of  the  emitter,  but  also  its 
linear  impedance,  and  also  in  a  certain  scale r  the  input  impedance  of 
the  element  in  the  grid  of  assigned  finite  dimensions  (with  uniform 
distribution  of  the  current  over  the  emitter) . 

As  shown  by  numerical  calculation  from  formula  (13)  ,  in  a  large 
grid  of  a  one-sided  emitter  with  current  distribution  diminishing 
towards  the  edges  of  the  type  (19)  ,  the  reactive  impedance  of  the 
element  is  comparatively  small.  At  the  same  time,  the  active  impedance 
of  the  element  in  this  grid  is  almost  constant  and  close  to  120  n . 

This  result  coincides  with  that  obtained  in  [3]  with  the  help  of 
the  concept  of  a  unified  cell.  It  indicates  the  basic  possibility  of 
broad-band  matching  of  elements  in  large  emitters.  We  note  that  by 
changing  the  c/d  ratio,  there  is  a  possibility  of  matching  the  average 
input  resistance  of  the  emitter  elements  with  feeder  of  any  assigned 
wave  impedance  W.  Thus,  for  example,  with  e2“^ *  w=?5  Ohm  and  one¬ 
sided  emission,  the  correlation  c/d=5  is  necessary. 

With  the  help  of  this  method  of  matching  the  elements  of  the 
emitter  through  selecting  the  distances  between  them,  in  certain  cases 
one  can  probably  obtain  wider-band  matching  than  with  the  help  of  the 
usually  employed  method  using  transforming  properties  of  the  elements 
themselves . 


Efficiency  of  the  Emitter  Placed  Near  the  Interface 

We  will  find  the  efficiency  of  the  emitter  that  is  placed  on  the 
interface  of  two  media.  An  approximate  expression  for  the  input 
impedance  of  a  large-sized  emitter  was  determined  by  formula  (17). 

With  the  help  of  (12)  ,  by  finding  the  impedance  of  the  radiation  under 
the  same  assumptions  ,  and  by  using  (8)  ,  we  obtair  the  efficiency  of 
the  large-sized  emitter 
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—LL 

Ve- + >r, 


(25) 


The  impedance  of  radiation  of  a  small  emitter  is  found  with  the 
same  assumptions  that  were  made  in  deriving  (18).  By  substituting 


in  (8)  ,  we  find  the  efficiency  of  the  small-sized  emitter 


1  = 


VJx 

**(>  ®i+l  **) 


(26) 


Figure  5  shows  by  solid  lines  the  dependences  of  efficiency  on 
disk  radius  for  media  with  different  parameters  that  were  obtained 
through  numerical  integration  of  accurate  expressions  (12)  and  (13). 

It  is  apparent  from  the  presented  curves  that  formula  (25)  yields 
satisfactory  accuracy  with  k2L  >  10.  On  the  other  hand,  formula  (26) 
yields  fairly  accurate  results  if  k2L  <  0.5. 

The  curves  in  fig.  5  determine  the  maximum  attainable  values  of 
efficiency  for  an  emitter  of  assigned  parameters  located  on  the  inter¬ 
face  if  medium  2  does  not  have  losses.  It  is  natural  that  when  there 
are  losses  in  medium  2,  the  efficiency  values  will  be  lower  than 
those  shown  in  fig.  5.  We  will  dwell  on  the  dependence  of  efficiency 
on  the  depth  of  submersion  of  the  emitter.  Figure  6  presents  the  de¬ 
pendences  of  efficiency  of  disk  emitters  of  varying  radius  on  the  depth 
of  submersion  for  a  medium  without  losses.  They  are  found  by  numerical 
integration  of  (12)  and  (13).  As  is  apparent  from  fig.  6,  for  large¬ 
sized  emitters  (k2L  >  6) ,a  drastic  peak  appears  in  the  dependence  of 
efficiency  on  depth.  Its  maximum  corresponds  to  the  equality  2k2h-ir. 
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For  small-sized  emitters,  this  peak  is  almost  not  pronounced. 

Figure  5  shows  through  the  dotted  lines  the  numerically  found  depen¬ 
dences  of  efficiency  with  optimal  deepening  on  the  disk  radius. 
Comparison  of  them  with  the  dependences  of  efficiency  with  zero  deepen¬ 
ing  shows  that  the  proper  selection  of  emitter  deepening  can  produce 
a  considerable  increase  in  the  efficiency  for  a  disk  of  sufficient 
radius  ^2^  <  6) .  It  can  be  shown  that  the  maximum  efficiency  which 
can  thus  be  obtained 

n  —  yi'~.  (27) 

1  +  Ke, 

Practical  realization  of  this  effect  is  difficult ,  first  of  all 
because,  as  is  apparent  from  fig.  6,  with  an  increase  in  permeability 
e 2»  the  width  of  the  peak  rapidly  contracts,  and  secondly,  because,  as 
shown  by  the  numerical  calculations  from  formulas  (12)  and  (13)  ,  the 
effect  disappears  when  the  imaginary  part  becomes  commensurate  with  the 
real  part. 

Conclusion 

The  averaged  input  impedance  of  square  (c/d=l)  elements  of  the 
emitter  which  are  connected  on  an  in  series-parallel  circuit  where  the 
emitter  radius  exceeds  two  wavelengths  in  a  low- conducting  medium,  is 
purely  active  and  equal  to  the  wave  impedance  of  the  medium  with 
one-sided  emission  and  half  of  the  wave  impedance  of  a  medium  with 
two-sided  emission. 

The  input  impedance  of  the  emitter  which  is  located  parallel  to 
the  interface  of  two  low-conducting  media,  with  a  change  in  the 
distance  from  it  to  the  interface ,  changes  in  limits  from 
the  quantity  of  wave  impedance  of  the  medium  with  great  permeability 
to  half  of  this  value. 

The  efficiency  of  an  emitter  located  on  the  boundary  of  two  low- 
conducting  media  whose  radius  exceeds  one  wavelength  in  the  medium 
with  lower  permeability,  is  approximately  proportional  to  the  root  of 
the  ratio  of  dielectric  permeabilities  of  the  media. 

The  efficiency  of  the  same  emitter ,  if  its  radius  is  less  than 
one-sixth  of  a  wave  in  the  medium  with  greater  permeability,  is 


approximately  proDortional  to  the  ratio  of  dielectric  permeabilities 
of  the  media  in  the  degree  3/2. 

With  deepening  of  the  emitters  of  sufficiently  large  radius,  one 
can  obtain  efficiency  in  the  low-conducting  medium  which  is  almost 
equal  to  a  unit. 

In  conclusion,  I  consider  it  my  duty  to  thank  Professor  B.  V. 
Braude  and  Professor  G.  A.  Lavrov  for  valuable  criticism  which  stimu¬ 
lated  the  publication  of  this  article. 


Appendix 

The  system  of  orthogonal  functions  (5)  can  be  the  solution  to  the 
vector  wave  equation  in  cylindrical  coordinates,  if  by  r , <p ,  and  z 
we  mean  the  coordinates  of  the  observation  point.  The  projections  of 
the  vectors  of  the  electric  and  magnetic  fields  on  the  axis  in 
the  region  between  the  interface  and  the  plane  passing  above  the 
sources  can  be  written  with  the  help  of  these  functions  as  follows: 

(!=  —  <»  0 


The  projections  on  the  r  axis  are  written  analogously  by  replacing 
the q?  signs  by  r.  The  flux  of  Pointing's  vector  passing  upwards 
through  the  plane  z=z^ , 

Pl  (2) (//|*  E'~H^  u rdrd  (29) 

By  substituting  (28)  and  by  using  the  conditions  for  orthogonality, 
we  obtain 
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as  well  as  the  obvious  correlations: 
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after  simple  transforms,  we  obtain  (1).  The  field  in  the  region  that 
lies  below  the  plane  which  places  below  the  sources  is  written  as: 


5  $l(frn+«>ti)Arn+(frn+«*ftn)A*)^di, 

/|*=— 00  0  * 

•o  o» 

V-yJ]  J((/2-„  +  «»fe)^  +  (fe  +  a3^)^l^-d/.  (32, 

n=-» o  * 

The  projections  on  the  r  axis  are  written  analogously  through  replace¬ 
ment  in  (32)  of  the  signs  ^  by  r.  By  substituting  (32)  into  the 
expression  for  the  flux  of  Pointing's  vector  which  passes  downwards 
through  the  plane  z=Z2 

~J  J  ("'rEU-H^^rdrd'p  (33) 

0  0 

and  by  using  (3)  and  (31),  we  obtain  (2). 
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ANALYSIS  OP  THE  PROPERTIES  OF  CONNECTED 
VIBRATOR  STRUCTURES 


N.  A.  Bey  and  N.  M.  Yakusheva 


Summary 


The  task  of  determining  the  phase  velocities  of  the  main  waves 
propagating  in  connected  infinite  vibrator  structures  is  solved 
by  the  method  of  integrated  equations.  The  obtained  correlations 
make  it  possible  to  link  the  phase  velocities  of  the  main  waves 
with  the  geometric  parameters  of  the  system. 


Introduction 


Small  transverse  dimensions',  small  connection  in  the  opposite 
direction,  and  the  possibility  of  forming  a  sector  beam  pattern  are  the 
advantages  of  antennas  of  axial  radiation,  and  in  particular,  director 
antennas  when  they  are  used  in  phased  antenn.a  arrays.  In  order  to 
solve  the  task  of  synthesizing  the  beam  pattern  of  the  emitter  for 
a  director- type  antenna  which  operates  in  an  array,  it  is  necessary  to 
know  the  size  of  the  linear  coupling  coefficient.  Its  determination 
is  reduced  to  calculation  of  the  phase  velocities  of  the  even  and 
odd  waves. 

It  should  be  noted  that  in  publication  [1]  which  also  covers  an 
analysis  of  the  properties  of  the  connected  vibrator  structures  ,  the 
obtained  correlations  are  very  cumbersome , and  only  the  properties  of 
the  structures  arranged  in  the  H-plane  with  equal  geometric  parameters 
are  numerically  analyzed. 

Derivations  of  the  Main  Correlations 

We  will  examine  a  general  case  of  connected  structures  which  are 
arranged  along  the  x-axis  and  are  shifted  along  the  y  and  z  axes  to 
distances  a  and  a  (fig.  1).  We  will  consider  the  vibrators  to  be 

y 

ideally  conducting  thin  tubes,  and  d/1  «  1  and  d/a  «  1.  The  boundary 

condition  will  be  the  equality  of  the  tangential  field  components  on 
the  vibrator  surface  to  zero.  As  the  conditions  for  infinit  y,  we  will 
use  the  condition  of  maximum  absorption.  Assume  that  all  the  vibrators 
are  excited  with  a  linear  change  in  phase  such  that  qax--phase  shift 
between  neighboring  vibrators  / _ ■ 

\  °x  «,/  ' 

Kn(z)  ~K*,  (;)pxp(i<7a(rt) 

Km  W  =  (z)  exp  (i  qa,m) , 

where  n  and  m--numbers  of  vibrators  of  first  and  second  channels. 

We  will  introduce  the  cylindrical  system  of  coordinates  z,  r,<|> 
with  beginning  in  the  center  of  the  zero  vibrator  of  the  first  channel. 
By  using  the  expression  obtained  in  [2]  for  the  vector  potential  near 
the  surface  of  the  zero  vibrator  of  the  first  channel ,  and  the  boundary 
condition,  after  certain  transforms  we  obtain  the  integrodif f erential 
equation  for  the  current  of  the  first  channel  vibrator 


-46- 


"nrrrrr' 


dV 

cfe? 


■  £  J  *X 


-H. 


u-/*)j  *;y;(«-) 

_ r'.  exp  [i «/(,,  -,,))  (rK- 


'  Vtx 


+ 


wk  c 

si :  £  j  <“x 


Wa— •  « — 0> 


x 


J  **  4  (  “P  (' 1  *  VT*  (i/I  -  a,  -  »;)  +  i  «  (?l  +<,,-*-)] 

F£  •  " 


Here: 


Figure  1. 
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U.%~  *—**“■  »  s~n  for  f^  and  m  for  f2* 

The  boundary  condition  has  the  appearance  y,  2(±J,.*)>=0. 

Analogously ,  an  equation  is  written  for  the  current  of  the  second 
channel  vibrator. 

In  [2],  the  integrodif ferential  equation  for  the  vibrator  currents 
is  solved  on  the  assumption  that  the  lengths  of  the  vibrators  differ 
little  from  the  resonance  (|k/i ,t-n  —  l«Sl).  Here  we  will  examine 


the  solution  for  the  case  where  detuning  of  the  vibrators  is  great,  i.e.  , 
'|*>/u  ~  |>|x!  In  a  similar  way  as  was  done  in  publication  [3], 


we  will  use  the  method  of  perturbations  ,  and  we  will  look  for  a  solution 
in  the  form  of  a  power  series  %'  xJ^+X*Jn+  •••  («i  --number 


of  channel  that  the  vibrator  belongs  to) . 

By  substituting  this  series  in  equation  (1)  and  the 

boundary  conditions  ,  we  obtain  a  system  of  differential  equations  with 
boundary  conditions.  The  solution  to  the  first  equation  will  be 
0. 

The  task  is  further  reduced  to  solving  the  system  of  equations: 


{zt)  K'J„  (z.)  =  ^  K,  (z).  Ju  (4  /,)  =  0 
J.',,  (Zj)"h  K*Jt  x(Zy)  ~  ~  j*  K  [GnqUl.  2)4-  ®*l«Vl2)l>  (i  U)—Q 
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J\7  (z.)  +  icV,,  (z.)  -  ^  Kt(z)-,  Ju  ( ±  U)  =  0 
J‘„(zJ  +  (2,)  «  {Jv  2)  4-  Gliq  (J, .  z)l,  Ju (4  /,)  =  0 

(3) 


The  solution  to  the  first  equations  in  systems  (2)  and  (3)  obtained 
in  publication  [3]  for  the  case  where  the  outside  emf  is  concentrated 
in  the  middle  of  the  vibrator  (z)  =  U6(z)\,  looks  like 


Ji  (*')  ^U'~^  lsin  (* 1  z'  ~tg  (Kl)  cos  (K  2,,l< 
We  will  determine  22q^l’  z^  an^  ^21  12^1*  : 
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The  solution  to  the  second  equations  of  systems  (2)  and  (3) 
sought  in  the  form 

(*!.:>)  Casin(irc)  * iJ4'1  £  j'  j,  . 


fl,  tfi  —  jo  -  -  X 


X - _ I  I ]  I  J _ **(*,,„  , 

0  /f)  J  -/i.»  [  cos J[  '  «p |  jj 

With  regard  for  the  boundary  conditions,  we  obtain: 
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After  transforms , (4)  acquires  the  appearance 
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v  T 1  a _ 1  cos  (*//,  n 
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At  the  point  z=0: 
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l  «pK^/-/,.2)  j 

The  solution  to  (2)  and  (3)  is  written  in  the  form 

A..-*',,., (6) 

After  substituting  (5)  into  (6) ,  we  obtain 


COS  (KflfO 


COS  (<//,  ,)  1*  f  _ COS  (n 

fi  f  m  i  1 114  eXPKM 

S  J  (t  -  <’) 


rO _ ' 


n,  nt~—o>  — 


In  a  similar  manner  as  was  done  in  [4],  we  find  the  equation  which 

links  the  phase  velocities  of  the  main  waves  which  are  propagated  in 

the  system  with  its  geometric  parameters  (with  accuracy  to  terms  on  the 
2 

order  of  x  ' 


r  COS  (*//,  ;)  j 

_tg(«/,j)+-?L-  y  f dti — 


x  14- 


cos  (m,I) 

mpW*  V 


In  particular,  for  a  single  structure,  (7)  acquires  the  appearance 

f.  «*(*"..?)  I* 


-tgK,)  |.?L  V  f  dti - MM  J  =0 


Results  of  Numerical  Calculation 

Figure  2  gives  the  relationship  (curve  1)  computed  according  to 
formula  (8),  and  for  comparison,  the  relationship  taken  from  [1]  (curve 
2)  which  also  characterizes  a  single  structure.  Figure  3  presents  the 
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curves  B=f(l/x)  of  even  and  odd  waves  which  propagate  in  the  connected 
structures  with  the  same  geometric  parameters,  arranged  in  the  H-plane. 
Figure  A  presents  analogous  curves  for  structures  in  the  E-plane.  Fig. 
5  presents  the  results  of  computing  the  relationship  B=f(l/A)  for  even 
(solid  line)  and  odd  (dotted  line)  waves  in  connected  structures,  and 
for  the  main  wave  of  a  single  structure  (dot-dash  line).  As  is 
apparent  from  fig.  5,  the  results  of  the  calculation  ffrom  the  formulas 
for  the  systems  made  of  slightly  detuned  vibrators  [2]  satisfactorily 
coincide  with  the  results  of  calculations  according  to  formulas  (7)  and 
(8)  only  for  1  >  (0. 215-0. 220) A . 
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ASSESSMENT  OF  AIRCRAFT  CAPACITIVE  CIRCUIT 
MEDIUM  WAVE  ANTENNAS  ACCORDING  TO  THEIR  EFFICIENCY 


E.  0.  Brudnyy  and  L.  Ya.  Il'nitskiy 


Summary 


This  work  presents  a  method  for  assessing  aircraft  capacitive 
circuit  medium  wave  antennas  according  to  their  efficiency. 

The  measure  of  efficiency  is  the  oroduct  of  the  antenna 
capacitance  times  the  active  altitude.  An  expression  is  obtained 
which  links  the  antenna  efficiency  to  its  geometric  dimensions. 

The  suggested  method  produces  the  optimal  dimensions  for  aircraft 
capacitive  circuit  antennas  direction  during  their  designing. 
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Capacitive  circuit  antennas. are  currently  used  as  nondirectional 
antennas  of  aircraft  medium-wave  range  radio  compasses.  However,  definite 
difficulties  arise  in  designing  these  antennas  in  the  selection  of 
the  optimal  geometric  dimensions  of  the  vertical  and  horizontal  circuit 
arms.  These  difficulties  are  due  to  the  lack  of  criteria  for  assessing 
the  operating  quality  of  the  circuit  antenna  associated  with  its 
geometric  dimensions. 

This  article  presents  a  method  for  assessing  the  capacitive 
circuit  antennas  according,  to  their  efficiency.  By  using  the  quasi¬ 
static  principle  to  analyze  the  operation  of  the  circuit  antenna  of 
medium  waves  [1-3],  one  can  draw  the  conclusion  that  the  main  measure 
of  its  efficiency  is  the  product  of  antenna  capacitance  times  the 
active  length.  It  is  quite  obvious  that  the  voltage  at  the  outlet  of 
the  radio  compass  receiver  will  be  proportional  to  this  product. 

The  antenna  efficiency  is  tied  to  the  geometric  dimensions  in  a 
quite  definite  manner. 

Aircraft  construction  usually  uses  a  circuit  antenna  which  con¬ 
sists  of  a  metal  tube  (or  conductor)  which  is  fastened  to  the  fuselage 
with  the  help  of  insulation  struts  and  closed  by  a  dielectric  cone. 

The  complete  capacitance  of  this  antenna  can  be  determined  from  formula 


[3] 


cA  = 


24.13  U  +  h) 


%  £[■+/  '-(£)’]! 


nf 


(I) 


where  l--length  of  horizontal  circuit  arm,m,  h--length  of  vertical 
circuit  arm,m  d--diameter  of  tube  (conductor)  of  which  the  circuit 
is  made ,  m. 

In  practice,  in  designing  the  aircraft  antennas  of  this  type, 
the  condition  2h  >>  d  is  always  fulfilled  .and  formula  (1)  can  be 
simplified 


24.13  U+h) 

•  Pf 


It  should  be  noted  that  the  real  capacitance  of  the  antenna  must 
somewhat  exceed  that  calculated  in  formula  (2) ,  since  here  the  calcula- 
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tion  docs  not  consider  the  continuous  ihsulator  which  increases  the 
total  capacitance  of  the  circuit'. 

Because  the  examined  antenna  operates  in  the  range  of  medium  waves , 
one  can  consider  the  distribution  of  current  on  the  vertical  circuit 
arm  to  be  uniform.  However  ,  this  conclusion  is  only  correct  with  a 
definite  correlation  between  the  lengths  of  the  horizontal  and  vertical 
antenna  arms . 

We  will  designate  through  hflQ  the  active  height  of  the  circuit 
antenna  which  is  installed  above  the  flat  counterweight,  then  [6] 

,  m 

b 

where  ^--current  in  the  antenna  base,  I (z) --distribution  function 
of  the  current  along  the  antenna. 

Assuming  that  the  horizontal  and  vertical  antenna  arms  are  made 
of  wire  of  equal  diameter,  and  that  the  following  condition  is  ful¬ 
filled  in  the  entire  working  range  of  the  antenna 

Xp«g  >  2X„ 

where  xpag --length  of  working  wave  of  antenna,  Xg--natural  length  of 
antenna  wave;  expression  (3)  can  be  written  in  the  form: 


where  W^-'Wave  impedance  of  vertical  circuit  arm, 

W^--wave  impedance  of  horizontal  circuit  arm,  k=^-. 

It  follows  from  an  analysis  of  formula  (4)  that  when  the  condition 
is  fulfilled 

with  accuracy  sufficient  for  practice,  one  can  consider  that 

(5) 

For  design  considerations,  the  height  of  the  aircraft  circuit 
antenna  is  usually  limited  by  the  condition: 

k <0,15  m. 
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then  the  equality  (5)  will  be  correct  for  the  antennas  where 

0,9-rl.O  m. 

Taking  into  consideration  that  expression  (5)  is  correct  only  for 
an  antenna  installed  above  the  flat  counterweight ,  we  will  introduce 
into  the  examination  h^-the  active  height  of  the  circuit  antenna  with 
regard  for  the  effect  of  the  aircraft  fuselage  ,  which  is  approximated 
by  the  surface  of  an  infinite  round  cylinder. 

It  is  known  from  publication  [4]  that 

hA  KyhAt,  Ms  (6) 


where  ky--ditnensionless  coefficient  of  antenna  elongation. 

For  the  circuit  antennas  ,  k^  depends  on  1 ,  increasing  proportionally 
to  the  length  of  the  horizontal  antenna  arm.  Thus,  for  example, 
for  an  antenna  with  horizontal  arm  1=1  m,  the  coefficient  of  elongation 
will  equal  2,  for  an  antenna  with  horizontal  arm  1=2  m,  the  coefficient 
of  elongation  will  equal  2.5,  and  for  an  antenna  with  1=3  m,  k^=3  [3], 


It  follows  from  expressions  (2)  and  (6)  that 

Sm*jSits±s  ,pf.» 


(7) 


where  S--efficiency  of  circuit  antenna. 


Figures  1-4  present  graphs  for  the  efficiency  of  aircraft  circuit 
antennas  for  different  values  of  their  design  dimensions.  The 
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relationships  (with  h  <<  A)  expressed  by  straight  lines  1,  are  computed 
for  the  circuits  with  length  of  the  horizontal  arm  1=3  m,  straight 
lines  2,  for  circuits  with  1=2  m,  and  straight  lines  3,  for  circuits 
with  1=1  m.  The  graph  calculations  were  made  according  to  formula  (7). 

If  the  minimum  value  of  efficiency  of  the  circuit  antenna  is 
assigned  based  on  the  conditions  to  guarantee  normal  operation  of 
the  radio  compass  (or  other  reception  device  that  the  antenna  must 
operate  with) ,  then  from  the  graphs  of  figs.  1-4  one  can  easily  deter¬ 
mine  the  geometric  dimensions  of  the  circuit  which  are  necessary  for  this. 
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The  suggested  method  was  experimentally  verified  [5].  The  experi¬ 
ment  was  conducted  with  a  capacitive  circuit  antenna  installed  on  the 
AN-24  aircraft,  and  had  geometric  dimensions  (h=0.01  m;  1=3.6  m,  d= 
0.007).  The  measurements  of  the  electric  parameters  Df  this  antenna 
yielded  the  following  values  (ho^0.2S  m  and  CA^76  pf ,  consequently, 
S«21.3  pf  x  m) .  The  minimum  efficiency  of  the  antenna  that  was 
necessary  for  the  normal  operation  of  the  radio  compass  was  selected 
as  equal  to  6  pf  m.  Then,  from  the  graph  1  in  fig.  3,  it  was  deter¬ 
mined  that  the  antenna  with  parameters  h=0.01  m,  1=2  m  and  d=0.007  m 
must  guarantee  the  assigned  efficiency.  After  shortening  the  hori¬ 
zontal  arm  of  the  antenna  to  2  m,  the  measurements  of  the  electric 
parameters  of  the  antenna  yielded  the  following  results:  A*  cti  0.23  m 
and  30  pf. 

Thus,  the  real  efficiency  of  the  antenna  became  equal  to  6.9  pf  m, 
while  that  obtained  from  the  graph  in  fig.  3,  7  pf  m.  Decrease  in 
antenna  efficiency  did  not  affect  the  operating  quality  of  the  radio 
compass  ,  while  the  weight  and  overall  dimensions  of  the  circuit  were 
considerably  reduced  [6]  , 

Taking  into  account  that  the  suggested  method  makes  it  possible 
to  obtain  the  optimal  dimensions  of  the  aircraft  capacitive  medium 
wave  circuit  antennas  directly  during  their  designing,  one  can  hope  that 
it  will  be  widely  employed. 
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